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Abstract. Magnetars are a subclass of neutron stars whose intense soft-gamma-ray bursts and quiescent X-ray 
emission are believed to be powered by the decay of a strong internal magnetic field. We reanalyze neutrino 
emission in such stars in the plausibly relevant regime in which the Landau band spacing AE of both protons and 
electrons is much larger than kT (where k is the Boltzmann constant and T is the temperature), but still much 
smaller than the Fermi energies. Focusing on the direct Urea process, we find that the emissivity oscillates as a 
function of density or magnetic field, peaking when the Fermi level of the protons or electrons lies about ~ 3kT 
above the bottom of any of their Landau bands. The oscillation amplitude is comparable to the average emissivity 
when AE is roughly the geometric mean of kT and the Fermi energy (excluding mass), i. e., at fields much weaker 
than required to confine all particles to the lowest Landau band. Since the density and magnetic field strength 
vary continuously inside the neutron star, there will be alternating surfaces of high and low emissivity. Globally, 
■ these oscillations tend to average out, making it unclear whether there will be any observable effects. 
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1. Introduction 



^ In recent years, evidence has been accumulating that soft gamma-ray repeaters and anomalous X-ray pulsars are mag- 
. . netars, strongly magnetized neutron stars powered by the decay of their magnetic field (see Woods & Thompson 2004 
^ ■ for a recent review) . The simplest estimate of their magnetic field strength relies on the magnetic dipole braking model 
k> \ usually applied to radio pulsars, which for these objects yields a very strong dipole field, ~ 10 14-15 G (Kouveliotou et 
■ al. 1998; Hurley et al. 1999). Contrary to the case of pulsars, their time-averaged luminosity (in the form of quiescent 
5^ , X-rays and soft gamma-ray bursts) by far exceeds their rotational energy loss, demanding another source of energy, 
which might be the decay of their internal magnetic field (Duncan & Thompson 1992; Paczynski 1992; Thompson & 
Duncan 1995, 1996), if this is substantially stronger than the external dipole field. 

The observable properties of these objects depend, among other things, on their field strength and their tempera- 
ture, whose evolution is coupled through the energy dissipated by the field decay and through the thermal modulation 
of the decay rates (Thompson & Duncan 1996). The most effective cooling mechanism for hot neutron stars is neutrino 
emission (see, e.g., Yakovlev & Pethick 2004), which also may to some extent control the magnetic field decay (Pethick 
1992; Goldreich & Reisenegger 1992). Particularly effective, if allowed, are the direct Urea processes, n — ► p + er + V e 
and p + e~ — ► n + i/ e , involving neutrons (n), protons (p), electrons (e - ), electron neutrinos (v e ), and electron an- 
tineutrinos (F e ). However, the matter inside a neutron star is highly degenerate. In the absence of an extremely strong 
magnetic field (far stronger than those discussed above; see Baiko & Yakovlev 1999), Pauli blocking, together with 
energy and momentum conservation, requires the triangle inequalities 

\PFp-PFe\<PFn<PFp+PFe (1) 

to be at least approximately satisfied for the direct Urea processes to be possible (pfi is the Fermi momentum of 
particle i, calculated in the non-magnetic approximation). Still, large uncertainties in the composition and equation of 
state of matter inside neutron stars cores prevent us from deciding whether this condition is fulfilled (e. g., Yakovlev & 
Pethick 2004). In this paper, we assume that it is, and study the effects of a magnetar-strength field on this simplest 



2 



M. Riquelme et al.: Neutrino emission in magnetized neutron stars 



neutrino emission process by accounting explicitly for the Landau-band structure of the proton and electron energy 
levels. We focus on the direct neutron /3 decay (n — > p + e~ + V e ) in chemical equilibrium. Other reactions should 
be affected in a qualitatively similar way, including modified Urea processes in chemical equilibrium or any reactions 
away from equilibrium, responsible for bulk viscosity (e. g., Finzi 1965; Haensel 1992; Reisenegger 1995). 

The magnetic field can make the emissivity change by two effects (Baiko & Yakovlev 1999). The first relates to 
the non-trivial spatial form of the wave functions of the charged particles. The second is an increased contribution 
to the emissivity from the highest Landau band occupied by electrons and/or protons. If the density or the magnetic 
field are changed, the emissivity oscillates (with an amplitude controlled by the temperature) as new Landau bands 
become populated or depopulated. The purpose of this work is to perform a detailed study of the second of these 
effects, quantifying its importance as a function of the different physical parameters involved. 

In order to justify the assumptions made in our calculations and to specify the situations in which this work applies, 
in §|2 we characterize the physical conditions we assume inside the star. In §[3] we describe the calculations performed 
and the results obtained. A general discussion of the relevance of the results is given in §0] whereas §0 lists our main 
conclusions. 



2. Physical Conditions 

It is known that a magnetic field B modifies the energy eigenstates and eigenvalues of electrons, protons and neutrons, 
whose masses we denote as m e , m p , and to„, respectively. The energies of free electrons and protons are given by 



E e = y ?7igC 4 + Pe Z c 2 + 2Behcn and E p = J rripC 4 + p pz c 2 + 2Behcm, (2) 

where c is the speed of light, e is the charge of the proton, h is the (reduced) Planck constant, and 

rn = I + - — ^~3pi where 1 = 0,1,2,... and a p = ±1. (3) 

The integer n and non-integer m label the Landau bands of electrons 1 and protons, respectively, a p is the doubled 
proton spin component along the direction of the magnetic field, g p = 2.79 is the gyromagnetic ratio of the proton, 
and pi Z is the component of the linear momentum of the particle i along the magnetic field direction. 

We note that, within each Landau band, the energy increases continuously with the absolute value of the longitu- 
dinal momentum component, \p% x \, starting from a bottom value corresponding to p lz — 0, which, for simplicity, we 
will refer to as the Landau level corresponding to this band. As the density increases in degenerate matter at a given 
magnetic field strength, new Landau bands get populated as the Fermi energy, Ept, increases, leading to oscillating 
thermodynamic functions (see, e. g., Dib & Espinosa 2001). Here, our main focus will be on similar oscillations in the 
neutrino emissivity. 

We ignore the effect of inter-particle forces on the energy levels of the charged particles. This is a good approximation 
for the electrons, which do not interact strongly. However, strong interactions do modify the eigenvalues of energy 
of the protons, whose Landau levels tend to broaden. We discuss below the importance of this broadening for the 
contribution to the emissivity from the highest Landau band occupied by protons. 

If we neglect strong interactions for the neutrons as well, the second triangle inequality in eq. would never be 
satisfied, and direct Urea reactions would be forbidden (e.g., Shapiro & Teukolsky 1983). Thus, for consistency, we 
allow for an arbitrary dispersion relation E n (p n ) for the neutrons, which takes their strong interactions into account 
and could in principle allow direct Urea processes to take place, even in the absence of a magnetic field. This relation 
defines an effective mass for the neutron, to* = PFn/[dE n /dp n ] PFn , which is generally somewhat smaller than the 
"bare" mass, m n . Here and throughout this paper, we ignore the relatively small dependence of the neutron dispersion 
relation on the spin projection along the magnetic field. 

In order to specify the physical conditions to be considered, we define the dimensionless quantities 

um = — ^ — , tom = - — , and r = — , (4) 

2Behc ' 2Behc ' 2Beh" V ' 

where the integer part of um corresponds to the highest Landau band populated by electrons in the zero-temperature 
limit. The same role is played for the protons by the largest value of m < ttlm compatible with eq. (0. We note that 
each of these can quantitatively be written as ~ 150p37 3 /-Bi6, where the number density of the relevant particles is 



Note that, except for the n = band, all the Landau bands of electrons are spin-degenerate, when QED corrections to 
i e = 2 are neglected. 
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P = P37 x 10 37 cm 3 and the magnetic field is B — Biq x 10 16 G. It is also useful to identify the energy difference 
between adjacent Landau levels for the particle i= e,p, 

AE { ~ ™*L, (5) 

1 /3 

which takes the values AE e ss 0.5Biq/ p^ 7 MeV for the (highly relativistic) electrons and AE P ss 0.06i?i6 MeV for 
the (approximately non-relativistic) protons. We assume that the electrons, protons, and neutrons are in chemical 
equilibrium and consider, for the calculation of the transition probabilities of the direct f3 decay, that the protons and 
neutrons are non-relativistic particles. 

Most importantly, we will see below that, in order for the Landau band structure to be noticeable, the magnetic field 
must be at least strong enough to make AEi ~ lOfcT, where k is the Boltzmann constant and T is the temperature. 
However, magnetic fields expected to exist in magnetars are still so weak as to make r, um, and tum 3> 1, which means 
that the electrons and protons populate many Landau bands. Writing T — Tg x 10 8 K, these conditions constrain the 
magnetic field to the range 

T 8 « B 16 « IOVst 3 - (6) 

While the upper bound is easily satisfied by all observational evidence, the lower bound is not obviously true in any 
known case. Indeed, the lower bound (here given for protons, although a less stringent condition for electrons also 

1/3 

exists: 0.1Tgp 3 y <C B\q) might not be satisfied even for magnetars, despite their very strong B fields, because the 
dissipation of their huge magnetic energy tends to keep their interiors fairly hot. If that is the case, the magnetic field 
effects we are studying here are negligible. Nevertheless, given the possibility that this regime may be applicable, we 
proceed to explore its consequences. 

3. Emission rate calculations 

We calculate the neutrino emissivity due to the direct (3 decay using the Weinberg-Salam-Glashow theory of weak 
interactions. We confirm Baiko & Yakovlev's (1999) eqs. (5) and (8) for the power radiated per unit volume, 

Q-- ^ 

where 



J d z p n d z p- v dp ez d Vpz E- v j n (\- f p )(l- f e )S(AE)6(A Pz )M, (7) 

n,l,(T v ,a n 



M = G% {2<^ [S^S^F^iu) + <^,-i<W^-i»] + [(! + 9AO P ) 2 F 2 n _ u (u) + (1 - g A a p ) 2 F^ »] } 



(8) 



u = c{p 2 n - (p ez +P P z) 2 }/ (2Beh), AE = E n — E p - E e - E„, Ap z = p nz - p ez - p pz , and fa is the Fermi-Dirac factor 
corresponding to the particle i. The Fermi weak coupling constant is Gf = 1.399 x 10 -49 erg cm 3 , and gA = 1.261 is 
the nucleon axial- vector coupling. The Laguerre functions F n j(u) are defined in the Appendix (and are taken to vanish 
when one of their indices is negative). The physical conditions discussed in 321 have already been used to neglect the 
antineutrino momentum, pp ~ kT/c <C PFm both in the definition of the variable u and in the momentum-conserving 
delta function. In principle, it contributes to the "thermal smoothing" discussed below, but to a lesser degree than 
the antineutrino energy, which we keep in our calculation. 

The integrals over directions of motion for antincutrinos and neutrons are easily done (the latter eliminates the 
momentum-conserving delta function), and the energies of antineutrinos, electrons, and protons can be used as new 
integration variables, yielding 

4Be r Pn EiE eEp f n {i- f p ){i- f^mi^M 



Q>= (0 ~ 96 ^ dE„ d Pn dE e dE p ^ (9) 

( 27r ) « c „,i,t~a n ,-y Jo Jo Ja ° Ja p J El - m 2 p c A - 2Behcm^/E 2 - m 2 e c A - 2Behcn 



where a e = ^m^c 4 + 2Behcn and a p = ^Jm^c 4 + 2BeTicm, and Q(u) is the usual Heaviside step function. The values 
of the new variable 7 = ±1 label the cases in which p ez and p pz have equal and opposite signs, respectively, affecting 

the value of u = p 2 n c 2 - (yW 2 - m 2 e c 4 - 2Behcn + jJe* - m 2 v c 4 - Wehcmj ^ /(2Behc). 

If the particles are in chemical equilibrium, the Fermi-Dirac factors suppress those reactions in which the difference 
between the energy of the reacting particles and their corresponding Fermi energies is much larger than kT . This fact 
allows to also change the integration variable for the neutrons, as p n dp n — m^dE n , and to set Ei = EFi in those terms 
where a variation of the energy of order kT makes very little difference. 
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Fig. 1. Magnetic-field dependent emissivity (normalized by the zero-field value) at low temperatures, as approximated 
by eq. (|llfl . for r = 100/& and um = tnM — 40/6, corresponding to constant particle number densities, p e — p p — 

0.25p„, and a magnetic field B = 6.8x 10 16 6 (p„/10 
effects at temperatures T ~ 10 9 K would strongly affect the emissivity. 



38 cm 3 ) 2 ^ 3 G. The shading identifies those regions where thermal 



For instance, in the Appendix it is shown that there is a range in which F mtn (u) behaves as an oscillat- 
ing function and that, outside this range, this function decays exponentially. The average oscillation period is ~ 
4-y/maxjn, m}/min{n, m} > 1, much larger than the variation in u due to an energy change ~ kT, Su ~ kT/ AE p -C 1. 

Therefore, in the variable u, we can safely set Ei = Epi, and obtain u = r — (y/riM — n + jy/rriM — fn) , where r, tim 
and rriM were defined in (4). 



With the factors \J E\ — m\ — 2Behcn and w E£ — — 2Behcm in the denominator of eq. @, we cannot always 

do the same. Since they could eventually vanish for the highest Landau band of each particle, a variation of the energy 
of order kT could become crucial in the emission rate we want to calculate. 

Let us first consider a situation in which this is not an issue, i.e., none of the Landau levels of the particles is too 
close to the respective Fermi energy, 



E Fe - \fr 



> kT and 



Ep p — \l rripC 4 + 2Behcm 



> kT (10) 



for all n and m. In this case, we can set E e = Ep e and E p — Ep p , obtaining 

q _ 4:57irm^m p EF e (kT) 6 \ - ^ - ^ - y^ 1 Q(u)M 

P ~ 80640?i 10 c 5 2^ Z- < 2L , 2-, ^j nM _ ny /m M - m ' 

7=±1 c„=±l CT P =±1 n=0 1=0 v v 

where the symbol [ J denotes the integer part, and l max = [m,M — \ + -^-\ . 

In Fig. we show a numerical calculation of this double sum for fixed densities of all particles, and a variable 
magnetic field. The true emissivity will deviate from this result (and the formal divergences will be smoothed, as 
discussed below) in those regions where the condition l|10l) does not hold. As an example, the shading in our graph 
identifies these regions for a temperature ~ 10 9 K. Still, it is instructive to note that the peaks come in "triplets", 
in which the central (and strongest) peak corresponds to a spin-degenerate electron Landau level coinciding with 
its Fermi energy, whereas the two flanking peaks correspond to similar coincidences for the (non-degenerate) proton 
Landau levels, for opposite spins. Their symmetric spacing arises automatically when the proton and electron densities 
are equal, as is the case in charge neutrality, when protons and electrons are the only charge carriers, but does not 
necessarily occur when other charge carriers (muons, mesons, or hyperons) are also present in the system. 

Much less conspicuous are the discontinuities caused by the step function Q(u) that ensures momentum conservation 
along the magnetic field. Since, as u — > 0, the Laguerre functions behave as F n ,i ~ u'™ - '' (see eq. IA.lt . a true 
discontinuity in the emissivity (rather than in one of its derivatives) arises only for terms with n = I (in which the 
spatial wave functions of protons and electrons coincide). Some of these are barely visible in Fig. 2] at b w 2.94, 
3.06, 3.25, and 3.40. These will also tend to be washed out by thermal effects, so we doubt that they could have any 
astrophysical significance. 

An average value for Q p can be obtained by making some approximations for the behavior of the Laguerre functions 
Fn,l when n, I, and \n — l\ are all large. We take the average expression for the functions F^ n (u) in the interval 
where they oscillate (eq. jA.llj ). and n (u) — where they undergo an exponential decay (see Appendix). Then, 
approximating the sums over n and I as integrals and neglecting the spin dependence of the particles' energies, we find 



(12) 
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I (X) 




-5 5 10 15 20 25 30 

Fig. 2. The function I(x), evaluated numerically from eq. (|14|l . 



which agrees with the unmagnetized case (Lattimer et al. 1991). 

The approximations leading up to eq. (|llf) are exact in the limit of vanishing temperature, but clearly break down 
at finite temperature when one of the Landau levels approaches the Fermi energy of the respective particle species, so 
the condition (|10|l is not satisfied for this level. The contribution of the latter to the total emissivity becomes very large 
(formally divergent, in the approximation considered above) and is no longer well-approximated by the corresponding 
term in the expression (|llf) . In order to determine the importance of this effect, we must compute the contribution 
of this Landau band independently and, then, compare it to Q p as given by eq. (|12f> . which still gives a reasonable 
approximation to the contribution of all the other bands. 

We focus on the case in which one of the Landau levels of the electrons, labelled as n (either h = [um\ or 
h = \jim\ +1), is close to their Fermi energy, i. e., \Ep e ~ yjm\c 4 + 2Behcn\ ~ kT, but the same does not happen for 
the protons. The contribution of the electron Landau band n to the emissivity, obtained from the appropriate term in 
eq. 0, summed over all proton Landau bands, and using the same approximations for F?[(u) as in equation ljl2|) . is 



| 2 Behcm* n m p Ep e ^ \ 2kT 



- GAl 2wv (fer) V ^ - ■*■) JW - (13) 



where x = (Ep e — y/m^c 1 + 2Behcn)/(kT) and 

m=r r x ^ t+u -^ dudt . (M) 

Jo Jo U 2(e x - U + l)(e'+"- x - 1) 

In Figure 2, we can see that this function has a peak at x = 3.08, where its value is 53.8, decreasing as I(x) w 
45 77T 6 /(5040v / ^) for x ^> 3, and exponentially for x < 0. The smooth shape of this function (left-right reversed, since 
a; is a decreasing function of B) is the true shape of the divergent peaks of Fig. 1 at finite temperature. 



4. Discussion 

We compare the emissivity contribution of the highest significantly populated electron Landau band to the unmagne- 
tized emissivity (which approximates the contribution from all the other Landau bands) by computing the ratio 

Qy _ I{x) Behc I(x) AE e 

Q, ~ 123 y/kTE Fe {El e - mlc 4 ) ~ 123 

As the density is increased or the magnetic field is decreased, the emissivity oscillates as new Landau bands start 
populating. Assuming, for definiteness, that the control parameter is Ep e , and ignoring variations of Ep p for now, 
these variations are roughly periodic with a period AE e , on which asymmetric peaks occur (sharply rising as the Landau 
level starts populating, slowly decreasing afterwards), with a total width at half maximum ~ lOfcT. If AE e 3> lOfcT, 
the fractional amplitude of the oscillations is ~ 0AAE e /\/kTEp e , becoming comparable to the total emissivity when 
the level spacing is the geometric mean of kT and the Fermi energy or, equivalently, B ~ 10 17 v / 78P37G, probably 
higher than is reached in magnetars. 

In the case of the protons, the situation is different. Since they interact strongly, their dispersion relation is not 
exactly the one shown in eq. Strong interactions, mainly with the neutrons, broaden the Landau levels of the 
protons. If this broadening is of order or larger than the spacing between the levels, the effect of emissivity oscillations 
given by the filling of their consecutive Landau bands tends to disappear. On the other hand, even if the protons did 
not interact strongly, the importance of the emissivity contribution of their highest populated Landau band would be 
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subdominant in comparison to the one of the electrons. This can be seen directly from eq. |15|). taking into account two 
facts. First, proton Landau bands are not spin degenerate, so the expression for the ratio Qv/Qa is reduced by a factor 
1/2. Second, their larger mass implies (at the same density) a smaller Landau-level spacing, yielding a more stringent 
condition to see substantial oscillations, and a smaller amplitude for these. It also implies that the last expression in 
eq. (|15|l would not be directly applicable, but would have to be replaced by a similar one, in which the mass term 
mpC 2 is subtracted from the Fermi energy. 

In a neutron star in diffusive equilibrium, the redshifted total chemical potential (including electrostatic potential 
energy) for each particle species is uniform throughout the region where the respective species is present. Since the 
redshift factor is a smooth function of position, so will be the local (or "intrinsic") chemical potentials, i. e., the Fermi 
energies considered in this work except at phase boundaries, where the electrostatic potentials can be discontinuous. 
We stress that this does not imply that the particle densities are smooth functions of position, but will rather have 
oscillations superimposed on the smooth increase with increasing Fermi energy, as illustrated, e.g., in Fig. 2 of Dib & 
Espinosa (2001). 

More directly to the point of this paper, even if the magnetic field is uniform, the smooth spatial variation of the 
Fermi energies will make the emissivity an oscillating function of position. Shells of larger and smaller emissivity will 
alternate within the star. If the magnetic field strength is truly uniform, these shells will be spherical, constant-density 
surfaces, otherwise they can be arbitrarily convolved. 

The observational relevance of this is not clear. As said above, in each "period" AE e , the emissivity peak has a 
width ~ lOkT and a fractional height ~ OAAEe/y/kTEFe, so its fractional contribution to the total emissivity over 
the corresponding period is - (lOfcT/ AE e )(0AAE e / \JkTE Fe ) = A^/kT/E Fe . The de generacy condition kT <C E Fe is 
amply satisfied in all neutron stars more than a few seconds old, so the contribution of the peaks to the total emissivity 
is very small. 



5. Conclusions 

We have shown that, in the regime in which the Landau level separation is intermediate between kT and the Fermi 
energies (excluding mc 2 ), the direct Urea emissivity is strongly enhanced when the Fermi energy of electrons or protons 
lies ~ 3kT above the bottom of any of the Landau bands. Thus, it will oscillate as a function of density or field strength. 
The oscillation amplitude becomes comparable to the total emissivity when the level separation is near the geometric 
mean of the other two energy scales. The required conditions may be present in magnetars, but this is not guaranteed. 
Even if it were the case, the large number of emissivity oscillations expected inside a given star could average out to 
give nearly the same global cooling rate as in an unmagnetized star, so their observational consequences are far from 
obvious. 



Appendix A: Some properties of the Laguerre functions 

For easy reference, in this Appendix we summarize some useful properties of the Laguerre functions, F m>n (u). These 
results were obtained by Kaminker & Yakovlev (1980, 1981). 

The Laguerre functions with indices m,n = 0, 1, 2, ... are defined as 



, A _ f y/n\/m\(-l) a e 5 a 5LJ(u), for m > n ( . . 

tm ^ u >~\ (~l) a F n , m (u) otherwise 1 ' 

where a = m — n and L^(u) correspond to the Laguerre polynomials, 

" r(n + q + l) (-u) k 

L » w = Ajr(Ha + i) k\(n-k)r (A ' 2) 



Defining f{u) = -^i{u — «_)(u — U+), where u± — (m + n + 1) ± yj Amn + 2(m + n) + 2, approximate expressions for 
F m . n (u) for different ranges of u in the limit n, m, |n — m\ S> 1, can be given as follows: 



— For u_ -u> 

F m , n (u) = (-1)^ {Airu^Jiujy 1 e u VfM-<t>. (A . 3) 

- Forw- U+ » (^^t-) 3 , 

F m .n(u) = (-1)" 1 {4nu^Jtujy i e -V7W-0 ; (A.4) 
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where 

u^+u+ ( {y/\ u ~ u +\ - V\ u - u -\} \ s /uZuT I {y/u-\u-u+\ + 

<t> = In ^ } — + In ^ r >— . (A.5) 

4 I — it_ I 2 I u\u+ — U-\ 

— For \u — U—\ Cii-,!^-!!-, 

F ro , n («) = (-1)^ f — p-i -V Ai ({u. - «} I^^Vl ■ (A.6) 



U- (u + — U- 

— For \u — u+\ <C u+ — tt_, 

*-<"> " '-""{ m^ 4 -.-) }'* 1 ( { "-" +) {^}*) • (A ' 7) 

where Ai(x) corresponds to the Airy function of the first kind (see, for example, Lcbcdcv 1972). The asymptotic 
behavior of this function, 

^ x ^°° (a - 8) 

I 7T a |£C] 4 COS(||X| 2 — j) X — » — 00, 

shows us that the Laguerre functions oscillate between u_ and u + and decay exponentially outside this range. 

— For u — u_ S> ( — — — ) and u + -u> ( — — — ] we have 

Vu+— u— J ' yu+— u— J ' 

F mj „ = |7m-v/-/(w)} 2 cos (uy/-f(u) - (pj , (A. 9) 

where 

V = 7-^ sin" 1 ± - ^4 — sin -t 1 ) — — (m + n + 2a + a ) . A.10 

In this work we use the average expression for n (u) in this range, 

(FLn) = {(« ~ ~ • (A.ll) 
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